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ABSTRACT 

FUSE spectra of the white dwarf stars G191-B2B, GD 394, WD 2211 - 495 
and WD 2331 — 475 cover the absorption features out of the ground electronic 
states of N I, N II, N III, O I and Ar I in the far ultraviolet, providing new 
insights on the origin of the partial ionization of the Local Interstellar Medium 
(LISM), and for the case of G191-B2B, the interstellar cloud that immediately 
surrounds the solar system. Toward these targets the interstellar abundances of 
Ar I, and sometimes N I, are significantly below their cosmic abundances relative 
to H I. In the diffuse interstellar medium, these elements are not likely to be 
depleted onto dust grains. Generally, we expect that Ar should be more strongly 
ionized than H (and also O and N whose ionizations are coupled to that of H 
via charge exchange reactions) because the cross section for the photoionization 
of Ar I is very high. Our finding that Ar I/H I is low may help to explain the 
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surprisingly high ionization of He in the LISM found by other investigators. 
Our result favors the interpretation that the ionization of the local medium is 
maintained by a strong EUV flux from nearby stars and hot gases, rather than 
an incomplete recovery from a past, more highly ionized condition. 

Subject headings: ISM: abundances — ISM: clouds — stars: white dwarfs — 
ultraviolet: ISM 



1. Introduction 

Our solar system is traveling through a low-density {uh ~ 0.1 — 0.3 cm~^), warm 
(T « 7000K), partiafly ionized gas cloud called the Local Interstellar Cloud (LIC) 
(Lallement et al. 1996) with a diameter of about 4 pc (Redfield & Linsky 2000). Near the 
LIC are similar clouds that give rise to absorption features with different velocities (Cry 
1996; Lallement 1996), known collectively as the Local Interstellar Medium (LISM). From 
the available evidence, the physical conditions in the LISM clouds are only sightly different 
from those of the LIC (Linsky 1996; Ferlet 1999). 

Presently, we do not have a clear understanding of the physical processes that are 
responsible for high fractional ionization of helium in the LISM. From observations that 
(n(He I)/n(H I)) = 0.07 (Dupuis et al. 1995), along with the knowledge that He/H = 0.10 
in all forms, it is evident that He is usually slightly more ionized than H. Also, helium 
shows less variability in how strongly it is ionized from one region to the next (Wolff, 
Koester, & Lallement 1999). While the EUV radiation from white dwarf and other stars 
can explain the observed fractional ionization of H (Vallerga 1998), these sources do not 
produce enough photons with energies above 24.6 eV to explain the ionization of He. To 
resolve this problem, two main proposals have been offered. One is that the LISM is not in 
a steady-state condition; it is returning from a much more highly ionized state produced 
by an energetic event in the recent past {t <^ 10® yr), such as the flash from a nearby 
supernova (Reynolds 1986; Frisch & Slavin 1996) or its shock wave (Lyu & Bruhweiler 
1996). The other possibihty is that the ionization of He is maintained in a steady state 
by the diffuse EUV radiation arising from conductive interfaces between the cloud edges 
and the surrounding, hot medium at T ~ lO^K (Slavin 1989; Slavin & Frisch 1998), or, 
alternatively, by recombination radiation from highly ionized but cooled gases that might 
surround us (Breitschwerdt & Schmutzler 1994). 

Sofia & Jenkins (1998) have proposed that the apparent abundances of the neutral 
forms of N, O and Ar, elements that are generally very lightly depleted in the interstellar 
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medium (if at all), can help to umavel the mystery of the He ionization. UV absorption 
lines in the spectra of white dwarf stars present good opportunities for studying these 
abundances in the LISM, since most of the brighter objects are well removed from major 
gas clouds elsewhere and remain immersed within a "Local Bubble" (diam. ~ 200 pc) of 
very hot, low density material surrounding the nearby cloud complex (Cox & Reynolds 
1987; Sfeir et al. 1999). We report here on initial spectra of four white dwarf stars taken 
with the Far Ultraviolet Spectroscopic Explorer (FUSE) (Moos et al. 2000; Sahnow et al. 
2000) taken at a resolution of 20 — 25 km s^^ (§§11—0), and we interpret the findings in the 
light of the Sofia & Jenkins proposal to obtain a better understanding on how the LISM is 
ionized (§|). 



2. Observations 

The stars, their locations in the sky, estimates for A^(H I) arising from the LIC only, and 
the wavelength ranges observed by FUSE are summarized in Table |l|. All of the stars are 
at distances that range between approximately 50—80 pc. The values of A^(H I)lic/A^(H I) 
(cf. Tables ^ and ^ indicate that the line of sight to G191-B2B contains mostly material 
from the LIC, while those toward the remaining three stars are dominated by other clouds 
in the LISM. 

Our spectra cover the important absorption features from N, O and Ar in their neutral 
forms. In addition, for some stars we recorded two ionized forms of nitrogen, N II and 
N III. While the N III feature conceivably might arise from photospheric absorption, we 
see no evidence of stronger absorption from N III in an excited fine-structure level - in the 
absence of saturation this feature would be twice as strong as the one we observed if all of 
the N III absorption came from the star's photosphere. Two spectral segments for one of 
the stars, WD 2211 — 495, are shown in Fig. |I[ 
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Table 1. Target Stars, Locations and Wavelengths Ranges 



Star^ 


i 


h 


logA^(H I) 


A Coveraj 






(deg) 


(deg) 


for the lAQ}" 


(A) 




G191-B2B 


156.0 


+7.1 


18.25 


987- lOJ 


\2 


GD 394 


91.4 


1.1 


17.68 


987 - lOJ 


12 


WD 2211 -495 


345.8 


-52.6 


16.73 


905 - lOJ 


12 


WD 2331 -475 


334.9 


-64.8 


17.12 


905 - lOJ 


12 



^Ordered according to increasing foreground A^(H I) - see 
Table | 

'^Estimate for the contribution to the hydrogen column that 
arises from the Local Interstellar Cloud surrounding our solar 
system, according to Redfield & Linsky (2000). 




1000 1020 1040 1060 1080 
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Fig. 1. — Far ultraviolet spectra of WD 2211—495 recorded by FUSE. Important features 
for our study are identified according to the numbers listed in the last column of Table |[ 
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Table 2. FUSE Measurements of Absorption Lines 





A 




WD 2211 - 495 


WD 2331 - 475 






(A) 




\W\ (mk)] 


\W\ (mk)] 


in Fief ffl 


NI ... 


952.303 


0.532^^ 


16.9 ±2.4 




1 




965.041 


0.589 




14.6 ±4.2 






954.104 


0.809 


4.8 ± 1.5 








964.626 


0.959 




23.4 ±4.1 






953.415 


1.098 




34.6 ±4.5 






963.990 


1.154 




22.4 ±4.2 






953.655 


1.376 


24.5 ±2.1 


34.4 ±4.1 


2 




953.970 


1.521 


23.2 ±2.0 


47.4 ± 4.8 


3 


N n . . . 


915.612 


2.180 


89. ± 15.= 


59. ± 20.^^ 


4 


N IF . . 


916.012 


2.180^^ 


-2.7 ±2.9*^ 


-2.9 ±5.3^ 




N IF* . 


916.701 


2.182^^ 


3.6 ±2.8*^ 


-9.0 ±5.7"^ 




N III . . 


989.799^ 


2.085 


64.5 ±4.0 


70.9 ± 4.0 


5 


N IIF . 


991.577^ 


2.085*^ 


26.1 ±2.3 


7.8 ±4.4 




OF... 


950.884 


0.176 


13.4 ±2.3 


29.1 ±5.6 


6 




929.517 


0.329 


27.3 ± 3.0 




7 




1026.476^ 


0.434 


4.3 ±2.1 


-0.01 ±3.0 






976.448 


0.509 


7.5 ±2.3 


58.2 ±6.0 






936.630 


0.534 


32.5 ±2.4 


40.6 ±4.7 


8 




948.685 


0.778 


50.4 ± 3.0 


69.9 ±6.9 


9 




1039.230'' 


0.980 


67.5 ± 1.5 


93.8 ± 1.8 


10 




971.738 


1.128 


55.5 ± 3.6 


124.0 ± 13.2 


11 




988.773' 


1.737^^ 


162.5 ±3.0 


184.6 ±7.5 


12 
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Table 2 — Continued 



Species 


A 


log(/A)^ 


WD 2211 - 495 


WD 2331 - 475 


Ident. 




(A) 




mimA)] 


[Wx{mA)] 


in Fig. H 


Ar P. . . 


1048.218'" 


2.440 


15.1 ± 1.0 


22.5 ± 1.8 


13 


Ar II . . 


919.781 


0.911 




0.3 ±6.1 





'^From a privately distributed update by D. C. Morton to his compilation published earlier (Morton 1991). 
^Two or more lines blended; line strength is for all of them combined. 

■^The large uncertainty is caused by an uncertain continuum in the vicinity of the closely spaced higher 
Lyman series transitions of H I. 

'^In principle, our upper limits for the column densities of N II in the excited fine structure levels, N II* 
and N II**, may be compared to lower limits for the amount of unexcited N II to arrive at a limit for the 
representative electron density (McKenna et al. 1996). Toward WD 2211 - 495 and WD 2331 - 475 we 
conclude that if T = 7000K, n(e) < 0.8 and 2.0 cm~^, respectively, but these values are higher than our actual 
expectations for the LISM based on C II*/C II and Mg I/Mg II seen elsewhere (Lallement &: Ferlet 1997; Wood 
& Linsky 1997; Jenkins, Gry, & Dupin 2000). 

'^G191-B2B: 51.0 =b 4.1mA; GD 394: 15.5 it 10.9mA. This line could have interference from an absorption by 
Si II A989.873, but it is not likely to be strong. 

fG191-B2B: -2.4 ± 4.0mA; GD 394: 15.4 ± 10.0mA. 

^G191-B2B: 1.3 ± 2.7mA; not measured for GD 394. This line consistently gave lower than expected 
absorption, based on the outcomes from other lines and the relative /-values. Hence, the published /-value for 
this line may be too large. For this reason, we disregarded indications from this line. 

'^G191-B2B: 35.5 ± 1.8mA; GD 394: 51.9 ± 3.4mA. 
'G191-B2B: 88.4 ± 3.9mA; GD 394: 139.4 ± 11.0mA. 

JThe line of Ar I at 1066.66A was not measured because of possible confusion with photospheric features of 
Si IV at 1066.61, 1066.64 and 1066. 65A. 

'"G191-B2B: 4.1 ± 1.5mA; GD 394: 12.0 ± 2.9mA. 
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3. Results 

The first and second columns in Table ^ show the species and the wavelengths of their 
transitions that we measured within the broad spectral coverages for WD 2211 — 495 and 
WD 2331 — 475, with the lines ranked according to their strengths expressed in terms 
of log(/A) (column 3) in each group. The fourth and fifth columns show the equivalent 
widths (in mA) and their errors arising from noise and uncertainties in the continuum and 
background levels. Our noise estimate was derived from the sizes of flux deviations from 
the adopted continuum level and thus included the effect of counting statistical fluctuations 
and small variations in the response of the detector (fiat-field corrections were not made). 
Results for the more limited set of lines that we could measure for G191-B2B and GD 394 
are given in the endnotes of the table. 

White dwarf stars can have narrow, photospheric features that could interfere with or 
masquerade as interstellar lines. This problem is probably more severe for the stars GD 394 
and WD 2211 — 495, which are known to be metal rich (Wolff et al. 1998). For example, 
the 952.3 A feature of N 1 in WD 2211 — 495 is probably enhanced by a stellar feature, since 
the observed equivalent width is far higher than those of other transitions of comparable 
strength. 

Table ^ summarizes our best estimates for the column densities of N 1, O 1, and Ar 1, 
along with the values for A^(H 1) reported by other investigators. A velocity dispersion 
b = 7km s^^ was derived for the O 1 lines toward WD 2211 — 495. The equivalent widths 
for the N 1 lines for this star are consistent with this 6- value. Toward WD 2331 — 475 
the N 1 lines give a best fit to 6 = 6 km s~^, while the O 1 lines are more consistent with 
b = 10 km s~^. This behavior probably arises from strong contrasts in physical environments 
for two or more velocity subcomponents, leading to different ratios of O 1 to N 1 in each 
case. For cases where the lines could not be seen or were marginally detected, we list 
upper limits for the column densities computed for the measurements plus their 2a errors, 
assuming that the lines have no saturation. For such large wavelength interval 

(~ O.3A) had to be considered (thus increasing the errors), because our wavelength scale 
was not very accurate. 

Our ultimate objective was to judge the column densities in the hght of a simple 
expectation based on some meaningful cosmic abundance scale. The entries within 
parentheses in Table ^ show how strongly deficient the neutral forms of these elements are 
with respect to their expected abundances relative to hydrogen, using the abundances in B 
stars as a standard. 
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Table 3. Neutral Atoms: Logarithms of Column Densities (cm ^) and Deficiencies' 



Species 


G191-B2B 


GD 394 


WD 2211 - 495 


WD 2331 - 475 


HI.... 
N I.... 


18.36^ 
13.90^ 
(-0.26) 


18.65l°if 
13.85 ±0.15^ 
(-0.60) 


18.76^ 
14.02 ±0.15 
(-0.54) 


18.93^ 
14.61 ±0.15 
(-0.12) 


1.... 


14.84'^ 
(-0.19) 


14.94 ±0.20^ 
(-0.38) 


15.321q;3 
(-0.11) 


15.45 ±0.1 
(-0.15) 


Ar I . . . 


< 12.44 
(< -0.42) 


12.70 ±0.1 
(-0.45) 


12.83 ±0.1 
(-0.43) 


13.06 ±0.18 
(-0.37) 



^Values in parentheses indicate by what logarithmic amount the abundances relative to H I deviate 
from their respective values in B stars, for which N = 7.80 (Cunha & Lambert 1994), O = 8.67 (Cunha 
& Lambert 1992) and Ar = 6.50 (Holmgren et al. 1990; Keenan et al. 1990) on a scale where H = 
12.00. 

(Vidal-Madjar et al. 1998). 
^ (Barstow et al. 1996). 
(Wolff et al. 1998); no errors stated. 

*^From an observation of the N I multiplet at 1200 A with GHRS on HST, reported by Barstow et al. 
(1996). While an arbitrarily large value of h gives a best fit to the lines, we expect that h ~ 10 km s~^ 
is more reasonable (and still gives an acceptable fit). 

^Based on our observations of the O I lines, but using h — 10km s~^ (see note e above). 

% — 6 km s~^ was adopted, as indicated by N I. The result drops to 13.01 \l b — 10 km s~^ (from 
1). 
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4. Interpretation 

The deficiency of Ar I for all four stars is of order —0.4 dex. We adopt the argument 
presented by Sofia & Jenkins (1998) that this is not caused by Ar depleting onto the surfaces 
of dust grains, but rather that neutral Ar is more easily photoionized than H. If some 
photons with E > 13.6 eV can penetrate a region, the H can be partially ionized. Under 
these circumstances the Ar (LP. = 15.76 eV) should be much more fully ionized because its 
photoionization cross section is about ten times that of H over a broad range of energies. 

N I and O I also have photoionization cross sections that are larger than that of H I, but 
the ionization fractions of N and O are coupled to that of H via resonant charge exchange 
reactions. For O II, the rate coefficient for charge exchange with H I (~ lO^^cm^s^^) (Field 
& Steigman 1971) is much larger than the coefficient of recombination with free electrons 
(4 X lO^^^cm^s"^ for T = 7000K), so the coupling of the O and H ionization fractions is 
very strong, unless there is appreciable ionization of O to higher stages. N II has a charge 
exchange rate that is only about twice the recombination coefficient (Butler & Dalgarno 
1979), so in regions where rig ^ ?t-h, N can start to behave more like Ar and show a 
deficiency of its neutral form. 

To illustrate the expected deviations for O I/H I, N I/H I and Ar I/H I away from 
the condition where all three elements are fully neutral. Fig. |^ shows the results of our 
calculations of photoionization equilibria for various depths of shielding by neutral H and 
He within a cloud. The detailed equations and how they are solved are described by Sofia 
& Jenkins (1998). At any point within the cloud the radiation field caused by hot stars 
was assumed to be stronger or weaker than Vallerga's (1998) combined stellar fluxes by 
a factor exp{[9 x lO^W^^ _ ^(jj I)][o-(A)h i + 0.07cr(A)He i]}, where the a(A)'s are the 
photoionization cross sections. We supplemented this fleld with Slavin's (1989) calculated 
flux from a cloud's conductive interface with much hotter gases.[^ 

Turning back to the observations, the results shown in Table |^ indicate that N and 
O are deflcient, but for G191-B2B and WD 2331 — 475 not as strongly as Ar. A mild 
deflciency in N can be understood if the hydrogen is mostly ionized, as is expected near the 
boundary of a cloud (see Fig. |^). The importance of ionization in reducing N I is reinforced 
by our lower limits (assuming no saturations of the proflles) for [logA^(N II), logA^(N III)] 
= [unavail., 13.68], [13.86, 13.78], [13.68, 13.82] toward G191-B2B, WD 2211 - 495 and 
WD 2331 — 475. Here, we can state that at least 38, 56 and 22% of all the nitrogen is in an 



We attenuated the strong emission peak near 70 A caused by Fe IX, Fe X and Fe XI, recognizing that 
McCammon et al. (2000) failed to observe these lines at a flux level far below the predicted one. This 
probably happens because Fe is depleted onto dust grains. 
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0.2 




Absorbing Column (cm ') 

Fig. 2. — Top panel: Logarithms of the expected deficiencies of the neutral fractions 
/o = 'ri(neutral)/n(total) of N, O and Ar, compared to those of H, caused by photoionization 
arising from the radiation by hot stars and a hot-gas conduction front surrounding a LISM 
cloud. These fractions are plotted as a function of the shielding depths [scaled to N(R I)] 
within a cloud at a uniform pressure p/k — 1.5 x 10^cm~^K and temperature T — 7000K. 
Local densities of neutral hydrogen and electrons (in cm"^) and the ratio of neutral H to He 
are shown in the bottom panel. 
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ionized form toward the respective stars. The underabundances of O are harder to explain, 
although they are reduced in magnitude by 0.17 dex if one adopts as a standard the 0/H 
found in the ISM by Meyer, et al. (1998) (8.50 on a logarithmic scale with H = 12.00) 
rather than the value of 0/H in B stars. Also, in view of the possibly large uncertainties 
in A^(H I) (except for G191-B2B), it may be better to adopt O as the reference element 
instead of H when we investigate the behavior of N 1 and Ar 1. 

For GD 394 and WD 2211 —495, it is difficult to understand why the deficiencies of N I 
are not much weaker than the corresponding ones of Ar 1. One possible explanation could 
be that a conduction front is not the dominant source of high energy radiation, but instead 
the ionization of N and Ar might be strongly infiuenced by 40 — 80 eV photons arising from 
recombinations to He II, which could conceivably be the dominant form of diffuse EUV 
radiation (Breitschwerdt, 2000) if much of the material in the Local Bubble has cooled to 
T ^ lO^K but has not yet recombined (Breitschwerdt & Schmutzler 1994). Over much of 
this energy range, N I is more easily ionized than Ar I if recombinations with free electrons 
are more important than charge exchange with H I [see Fig. 3 of Sofia & Jenkins (1998)]. It 
will be interesting to see if detailed calculations confirm this conjecture. 

A flux of ionizing photons with E > 24.6 eV and sufficient intensity to create the 
steady-state ionization of He in the LISM has not yet been observed directly, but our 
observations showing a significant reduction in Ar I favor its existence. If, by contrast, we 
had found a normal amount of Ar I, we might have concluded that the gas was probably 
more highly ionized in the past and has not yet fully recombined. Here, a near equality 
of Ar I/H I to the expected cosmic ratio would have arisen from their virtually identical 
recombination coefficients. Ionizations to levels higher than the singly ionized forms are not 
important in this picture, since their recombination rates are significantly faster (oc Z"^). 

We thank T. M. Tripp and C. Howk for their review and comments on an early draft 
of this paper. 
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